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Abstract 

The present knowledge on the association of single nucleotide polymorphisms (SNPs) of lysyl oxidase-like 1 (LOXL1) with 
pseudoexfoliation syndrome (PEXS) and pseudoexfoliation glaucoma (PEXG) is controversial and inconclusive. This meta- 
analysis sought to derive a more precise estimation of the effects of LOXL1 SNP loci (rsl 048661, rs3825942, and rs21 65241) 
on PEXS/PEXG. Literature searches were conducted on the PubMed, EMBASE, ISI Web of Science, and Cochrane Library 
databases through October 2013. Twelve studies describing 1810 cases and 1790 controls met the inclusion criteria. The 
strengths of the associations found through the meta-analysis were assessed with pooled odds ratios and their 95% 
confidence intervals (CI). A meta-regression analysis was also used to examine the influence of the study and population 
characteristics. The results indicated that rsl 048661 TT carriers had 92.1% and 40.4% less risk of developing PEXS/PEXG than 
did the controls in the Caucasian and Asian populations, respectively. Carriers of rs3825942 AA or rs21 65241 CC also had 
significantly less PEXS/PEXG susceptibility than did the non-carriers. Meta-regression showed that in Caucasians, the male 
proportion (slope: 0.272; 95% CI: 0.167-0.376; P = 0.0001) and mean age (slope: 0.796; 95% CI: 0.375-1.217; P = 0.0002) of the 
PEXS/PEXG subjects correlated positively with the effect of rs3825942 on PEXS/PEXG susceptibility. The meta-analysis 
suggested that LOXL1 rsl 048661 TT, rs3825942 AA, and rs21 65241 CC were associated with a reduced risk of developing 
PEXS/PEXG. 
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Introduction 

Pseudoexfoliation syndrome (PEXS) is an age-related systemic 
disorder that is the most common cause of secondary glaucoma 
worldwide and the most frequent cause of unilateral glaucoma 
[1,2]. In addition, PEXS progresses to pseudoexfoliation glaucoma 
(PEXG), which responds poorly to medical therapy in comparison 
to other types of glaucoma and can lead to the rapid progression of 
optic nerve damage [3], 

Despite its worldwide distribution, there is a clear tendency for 
PEXS to cluster geographically and in certain racial or ethnic 
subgroups [4]. Furthermore, PEXS has a strong familial associ- 
ation [5] . The underlying causes of the different prevalence rates 
between age-matched geographical and ethnic populations remain 
unknown, but appear to be related to variation in genetic 
background [5,6]. Lysyl oxidase-like 1 (LOXL1) is a member of 
the lysyl oxidase gene family and is essential to the biogenesis of 
connective tissue [7] ; it encodes an extracellular copper-dependent 
amine oxidase that catalyzes the first step in the formation of 



crosslinks in collagens and elastin [7,8]. A highly conserved amino 
acid sequence at the C-terminus appears to be sufficient for amine 
oxidase activity, suggesting that all family members may retain this 
function. A fibrillar, proteinaceous substance, is produced in 
abnormally high concentrations within the ocular tissues of 
patients with PEXG, and LOXL1 may be relevant to its formation 
[9,10]. 

A genome-wide association study identified three common 
single nucleotide polymorphisms (SNPs) in the LOXL1 gene on 
chromosome 15q24. 1 that were strongly associated with pseu- 
doexfoliation syndrome [11]. The LOXL1 polymorphisms 
included one intronic SNP, rs2 165241, located within the first 
intron, and two non-synonymous coding SNPs, rsl 048661 and 
rs3825942, located within the first exon. The association of 
LOXL1 with PEXS/PEXG has recently been confirmed in 
several different populations [12,13,14,15,16]. However, the 
reported associations are controversial and inconclusive due to 
factors including the limited sample sizes,different ethnicities, and 
genotyping procedures [17,18]. Believing a meta-analysis would 
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Figure 1. Flowchart of the study selection. SNPs, single nucleotide polymorphisms. 
doi:10.1371/journal.pone.0090331.g001 



provide more credible evidence by systematically summarizing the 
existing data, we gathered eligible studies to investigate the 
association between the LOXL1 gene polymorphisms and 
susceptibility to PEXS/PEXG. 

Materials and Methods 

Search strategy 

The PubMed, EMBASE, ISI Web of Science, and the 
Cochrane Library databases were electronically searched for 
case-control studies published through October 2013 that 
examined the association of the LOXL1 gene polymorphisms 
with the PEXS/PEXG susceptibility. The search strategy was 
based on a combination of "(lysyl oxidase-like 1 OR LOXL1) 
AND (gene OR variants OR polymorphism OR alleles OR 
mutation) AND (pseudoexfoliation syndrome OR pseudoexfolia- 
tion glaucoma)". We also manually searched references in key 
articles. The language was limited to English. 

Selection criteria 

The inclusion criteria were as follows: (a) evaluation of the 
association of the LOXL1 gene polymorphisms with the PEXS/ 
PEXG risk; (b) case-control studies; (c) sufficient published data for 
estimating an odds ratio (OR) with a 95% confidence interval (CI); 
and (d) PEXS was diagnosed as the presence of pseudoexfoliative 
material on the anterior lens capsule after maximal pupil 
dilatation. PEXG was diagnosed if patients had typical features 
of PXFS and all of the following: an initial intraocular pressure of 
at least 22 mm Hg, glaucomatous optic disc changes, visual field 
defects consistentent with optic nerve damage, and no evidence of 
other conditions causing secondary glaucoma [19]. The study 
authors were contacted for supplemental data if the information 



was not available in the publication. The studies with overlapping 
patient samples were excluded, and only the studies with the larger 
numbers of patients were included. Two of the authors indepen- 
dently identified and reviewed each relevant paper, and disagree- 
ments were reconciled through group discussion. 

Data extraction 

To show the relationship between the LOXL1 gene polymor- 
phisms and PEXS/PEXG risk, the most strongly and indepen- 
dently associated SNPs were selected for the analysis (rs 1048661, 
rs3825942, and rs2 165241). Information regarding the following 
aspects was retrieved from each study according to a fixed 
protocol: study design; geographical location; population ethnicity, 
definition and numbers of cases and controls; DNA extraction and 
genotyping methods; frequency of the genotypes; mean age of the 
patients; and proportion of the patients who were male. When the 
studies included subjects of more than one ethnicity, the genotype 
data were extracted separately for each ethnic group. The genetic 
equilibrium of the LOXL1 gene for the control group of each 
study was evaluated by testing for Hardy-Weinberg equilibrium 
(HWE) using chi-square analyses [20] . A state of disequilibrium 
was defined as a P value <0.05. 

Statistical analysis 

A summary OR with a 95% CI was calculated to assess the 
strength of the association of the LOXL1 gene polymorphisms 
with the PEXS/PEXG risk. The OR of each study was first 
calculated in a 2x2 table. Pooled ORs for the risk were then 
calculated for the allele frequency comparison and the additive, 
dominant, and recessive models. The between-study heterogeneity 
was evaluated with a 0_ statistic, and a P value <0.1 was 
considered statistically significant [21]. If the P value was >0.1, a 
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Table 1. Main characteristics of all studies included in the meta-analysis. 





First author (year) 


Population 
ethnicity 


LOXL1 dpSNP rsID, Allele PEXS/PEXG Controls 


Case males 
n (%) 


Case mean age 
(year) 


HWE* 


Challa (2008) 


Caucasian 


rs 1048661 T/G 








YES 






rs3825942 A/G 


50 235 


39 (78.0) 


74.0 


YES 






rs2165241 C/T 








YES 


Fan (2008) 


Caucasian 


rs 1048661 T/G 








NO 






rs3825942 A/G 


199 116 


84 (42.2) 


75.0 


NO 






rs2165241 C/T 








YES 


Ozaki (2008) 


Asian 


rs 1048661 T/G 








YES 






rs3825942 A/G 


209 1 72 


67 (32.1) 


78.0 


YES 






rs2165241 C/T 








YES 


Pasutto (2008) 


Caucasian 


rsi 048661 T/G 








YES 






rs3825942 A/G 


726 412 


312 (43.0) 


77.1 


NO 






rs2 165241 C/T 








NO 


Ramprasad (2008) 


Asian 


rsi 048661 T/G 


52 97 


27 (51.9) 


68.9 


YES 






rs3825942 A/G 








YES 


Lee (2009) 


Asian 


rsi 048661 T/G 


62 171 


30 (48.4) 


74.7 


YES 






rs3825942 A/G 








YES 


Abu-Amero (2010) 


Asian 


rsi 048661 T/G 


93 101 


61 (65.6) 


72.3 


YES 






rs3825942 A/G 








YES 


Malukiewicz (201 1) 


Caucasian 


rs 1048661 T/G 








YES 






rs3825942 A/G 


36 30 


9 (25.0) 


73.0 


YES 






rs2165241 C/T 








YES 


Jaimes (2012) 


Latin American 


rs 1048661 T/G 








NO 






rs3825942 A/G 


102 97 


NA 


74.8 


YES 






rs2165241 C/T 








YES 


Micheal (2012) 


Asian 


rs 1048661 T/G 


128 180 


69 (53.9) 


47.3 


YES 






rs3825942 A/G 








YES 


Metaxaki (2013) 


Caucasian 


rsi 048661 T/G 








YES 






rs3825942 A/G 


48 52 


43 (49.4) 


77.5 


NO 






rs2165241 C/T 








YES 


Park (2013) 


Asian 


rsi 048661 T/G 








YES 






rs3825942 A/G 


110 127 


53 (47.3) 


71.6 


YES 






rs2165241 C/T 








YES 


HWE, Hardy-Weinberg equilibrium; LOXL1, lysyl oxidase-like 1; NA, not available; PEXG, pseudoexfoliation glaucoma; PEXS, pseudoexfoliation syndrome; SNP, single 



nucleotide polymorphism. 

* The genetic equilibrium of the LOXL1 gene for the control group of each study was evaluated by testing for HWE using chi-square analyses. Disequilibrium was 
defined as P<0.05. 

doi:1 0.1 371 /joumal.pone.0090331 .t001 



fixed-effect model was used for the meta-analysis; otherwise, a 
random-effect model was used. The significance of the pooled OR 
was determined with the Z-test, and P<0.05 was considered 
statistically significant. 

To consider the possible sources of heterogeneity, we stratified 
the studies by ethnicity and repeated the analysis separately for 
each group. Furthermore, we performed a meta-regression 
analysis to assess the influence of the population characteristics 
[22]. We also performed sensitivity analyses, serially excluding 
studies to determine the sources of heterogeneity and assess the 
stability of the results. The Begg funnel plot asymmetry 
wasassessed with Egger linear regression tests, a linear regression 
approach that measures funnel plot asymmetry on the natural 
logarithm scale of the OR. The significance of the intercept was 



determined by the /-test suggested by Egger, with P<0.05 
considered to represent statistically significant publication bias 
[23]. 

All calculations were performed using the Comprehensive 
Meta-Analysis computer program version 2 (Biostat, Englewood, 
NJ, USA). 

Results 

Characteristics of the articles in the analysis 

The first search retrieved 83 articles. After eliminating the 60 
studies not related to the topic, 23 potentially relevant studies were 
identified for further evaluation. In these studies, the titles or 
abstracts indicated that they evaluated the association between the 
LOXL1 gene polymorphisms and PEXS/PEXG susceptibility; 
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Table 2. Summary of pooled odds ratios for the association of rsl 048661 and PEXS/PEXG in the meta-analysis. 





Ethnicity 


Genetic model 




OR (95% CI) 


P 


Heterogeneity test 


Publication bias 














12 


P (Q-test) 


P value 




Total 


allele frequency 
comparison 


Tvs. G 


0.942 (0.617-1.436) 


0.780 


97.2 


0.001 


0.400 






additive model 


TG vs. GG 


0.545 (0.421-0.705) 


0.001 


71.6 


0.001 


0.185 






additive model 


TT vs. GG 


0.294 (0.149-0.579) 


0.001 


82.9 


0.001 


0.921 






dominant model 


TT+TG vs. GG 


0.538 (0.377-0.768) 


0.001 


83.7 


0.001 


0.959 






recessive model 


TT vs. TG+GG 


0.393 (0.196-0.791) 


0.009 


88.1 


0.001 


0.788 




Caucasian 


allele frequency 
comparison 


Tvs. G 


0.439 (0.154-1.251) 


0.123 


96.5 


0.001 


0.859 






additive model 


TG vs. GG 


0.678 (0.361-1.273) 


0.227 


67.5 


0.027 


0.146 






additive model 


TT vs. GG 


0.067 (0.014-0.327) 


0.001 


82.2 


0.001 


0.829 






dominant model 


TT+TG vs. GG 


0.339 (0.174-0.661) 


0.001 


80.3 


0.002 


0.910 






recessive model 


TT vs. TG+GG 


0.079 (0.012-0.513) 


0.008 


88.0 


0.001 


0.775 




Asian 


allele frequency 
comparison 


Tvs. G 


1.056 (0.262-4.261) 


0.939 


97.8 


0.001 


0.748 






additive model 


TG vs. GG 


0.386 (0.283-0.528) 


0.001 


0 


0.998 


0.439 






additive model 


TT vs. GG 


0.447 (0.250-0.800) 


0.007 


47.2 


0.099 


0.304 






dominant model 


1 1 + 1 U VS. UU 


n Ton in i n 1 n\ 


0.001 


75.8 


0.002 


0.286 






recessive model 


TT vs. TG+GG 


0.596 (0.228-1.559) 


0.291 


48.9 


0.063 


0.236 




Latin American 


allele frequency 
comparison 


Tvs. G 


1.096 (0.673-1.788) 


0.712 


0 


1.000 


# 






additive model 


TG vs. GG 


1.889 (0.989-3.609) 


0.054 


0 


1.000 


# 






additive model 


TT vs. GG 


0.480 (0.141-1.636) 


0.241 


0 


1.000 


# 






dominant model 


TT+TG vs. GG 


1.452 (0.804-2.619) 


0.216 


0 


1.000 


# 






recessive model 


TT vs. TG+GG 


0.399 (0.119-1.342) 


0.138 


0 


1.000 


# 




CI, confidence interval; I 2 , inconsistency index; OR, odds ratio; PEXG, pseudoexfoliation gl 


aucoma; PEXS, pseudoexfoliation syndrome; SNP, 


single nucleotide 





polymorphism; vs., versus. Bold text: significant odds ratio (P<0.05) and significant between-study heterogeneity (P<0.1). 
# Publication bias could not be tested because a minimum of 3 studies were required. 
doi:1 0.1 371/journal.pone.0090331 .t002 



however, some studies did not meet all of the study inclusion 
criteria. Ultimately, 1 1 studies were excluded because they were 
undesirable article types (review or letter; n = 8), had overlapping 
patient samples (n= 1) [24], insufficient data (n = 1) [25], or were 
not related to the three SNPs (n = 1) [26] . Finally, 1 2 studies with a 
total of 1810 cases and 1790 controls were included in the analysis. 
Figure 1 shows the flowchart for the study selection. 

Table 1 shows the primary characteristics of the studies included 
in this meta-analysis. All of the studies used healthy control 
subjects, and the sample sizes ranged from 66 to 1 133. The studies 
originated from Europe (n = 5) [12,14,27,28,29], Asia (n = 6) 
[13,15,17,18,30,31], and Latin America (n= 1) [16]. Two studies 
extracted only the allele frequencies [12,17]. The study by Park et 
al. [18] represented the data as (minor homozygosity + heterozy- 
gosity) /(major homozygosity). There were no significant differ- 
ences between the case and control subjects with respect to age 
distribution. 

Quantitative synthesis 

Association of the rsl048661 T/G polymorphism with 
PEXS/PEXG. Twelve studies evaluated the association be- 
tween the rsl 048661 T/G polymorphism and the risk of 
developing PEXS/PEXG. For the combined group data, signif- 
icant associations between rsl 048661 and susceptibility to PEXS/ 
PEXG were observed for the additive, dominant, and recessive 



models, but not for the allele frequency comparison (Table 2). In a 
subgroup analysis performed by ethnicity, the pooled OR for the 
Caucasian population indicated a significandy decreased risk for 
PEXS/PEXG in the dominant (Figure 2), recessive, and TT 
versus GG additive models (Table 2). For the Asian population, 
the dominant (Figure 2) and two additive models showed a 
significant association between rsl 048661 and a reduced suscep- 
tibility to PEXS/PEXG, while the allele frequency comparison 
and recessive model did not (Table 2). No significant associations 
were found in the Latin American population. 

Association of the rs3825942 A/G polymorphism with 
PEXS/PEXG. The twelve studies contained data regarding the 
association of the rs3825942 A/G polymorphism with the 
susceptibility to PEXS/PEXG. In comparison with the control 
group, the association between rs3825942 and decreased suscep- 
tibility to PEXS/PEXG was significant in all genetic models 
(Table 3). Furthermore, in the subgroup analysis by ethnicity, the 
genetic models showed a significant association between 
rs3825942 and a reduced susceptibility to PEXS/PEXG in the 
Caucasian and Asian populations (Table 3; Figure 3). In the Latin 
American population, the allele frequency comparison and the 
dominant model showed significant associations (Table 3) and a 
reduced risk (Figure 3). 

Association of the rs2165241 C/T polymorphism and 
PEXS/PEXG. Eight studies contained data for the rs216524 C/ 
T polymorphism. The association between rs2 16524 and suscep- 
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Model 


Group by 


Study name 




Statistics for each study 






Ethnicity 




UQUS 


Lower 


Upper 












ratio 


limit 


limit 


Z-Value 


p-value 




Asian 


Ramprasad, 2008 


0.468 


0.236 


0.928 


-2.172 


0.030 




Asian 


Lee, 2009 


0.429 


0.220 


0.834 


-2.493 


0.013 




Asian 


AbuWmero, 2010 


0.378 


0.202 


0.706 


-3.051 


0.002 




Asian 


Micheal,2012 


0.311 


0.192 


0.504 


•4.744 


0.000 




Asian 


ParK2013 


0.050 


0.019 


0.131 


-6.105 


0.000 


Fixed 


Asian 




0.314 


0.236 


0.417 


-7.969 


0000 


Random 


Asian 




0.280 


0.154 


0.510 


-4.160 


0.000 




Caucasian 


Fan, 2008 


0.137 


0.081 


0.232 


-7.415 


0.000 




Caucasian 


Pasutto, 2008 


0.351 


0.273 


0.450 


-8.221 


0.000 




Caucasian 


Malukewicz, 2011 


0.483 


0.157 


1.482 


-1.273 


0.203 




Caucasian 


Metaxak. 2013 


0.778 


0.334 


1.810 


-0.583 


0.560 


Fixed 


Caucasian 




0.319 


0.258 


0.396 


-10.454 


0.000 


Random 


Caucasian 




0.339 


0.174 


0.661 


-3.175 


0.001 




Latin American 


Jaimes.2012 


1.452 


0.804 


2.619 


1.237 


0.216 


Fixed 


Latin American 




1.452 


0.804 


2.619 


1.237 


0.216 


Random 


Latin American 




1.452 


0.804 


2.619 


1.237 


0.216 


Fixed 


Overall 




0.357 


0.303 


0.421 


-12.280 


0.000 


Random 


Overall 




0.538 


0.377 


0.768 


-3.417 


0.001 



Odds ratio and 95% a 




Figure 2. Forest plots for associations between rsl 048661 T/G (TT+TG vs. GG) and PEXS/PEXG stratified by ethnicity. The first author 
of the study and year of publication are shown for each citation. The squares and horizontal lines correspond to the study specific odds ratio (OR) and 
95% confidence interval (CI), respectively; the diamonds represent the pooled OR and 95% CI. 
doi:10.1371/journal.pone.0090331.g002 



Table 3. Summary of pooled odds ratios for the association of rs3825942 and PEXS/PEXG in the meta-analysis. 





Ethnicity 


Genetic model 




OR (95% CI) 


P 


Heterogeneity test 


Publication bias 












12 


P (Q-test) 


P value 


Total 


allele frequency 
comparison 


A vs. G 


0.153 (0.104-0.225) 


0.001 


86.6 


0.001 


0.628 




additive model 


AG vs. GG 


0.153 (0.096-0.244) 


0.001 


42.6 


0.083 


0.296 




additive model 


AA vs. GG 


0.101 (0.052-0.198) 


0.001 


0 


0.933 


0.643 




dominant model 


AA+AG vs. GG 


0.015 (0.006-0.038) 


0.001 


86.4 


0.001 


0.365 




recessive model 


AA vs. AG+GG 


0.128 (0.066-0.250) 


0.001 


0 


0.913 


0.653 


Caucasian 


allele frequency 
comparison 


A vs. G 


0.277 (0.075-1.031) 


0.056 


93.7 


0.001 


0.883 




additive model 


AG vs. GG 


0.164 (0.057-0.473) 


0.001 


75.5 


0.007 


0.767 




additive model 


AA vs. GG 


0.085 (0.039-0.189) 


0.001 


0 


0.775 


0.799 




dominant model 


AA+AG vs. GG 


0.030 (0.006-0.149) 


0.001 


91.5 


0.001 


0.971 




recessive model 


AA vs. AG+GG 


0.106 (0.048-0.233) 


0.001 


0 


0.777 


0.790 


Asian 


allele frequency 
comparison 


A vs. G 


0.148 (0.099-0.223) 


0.001 


0 


0.587 


0.520 




additive model 


AG vs. GG 


0.156 (0.092-0.264) 


0.001 


0 


0.885 


0.068 




additive model 


AA vs. GG 


0.139 (0.036-0.535) 


0.004 


0 


0.778 


0.378 




dominant model 


AA+AG vs. GG 


0.019 (0.006-0.064) 


0.001 


82.1 


0.001 


0.185 




recessive model 


AA vs. AG+GG 


0.192 (0.050-0.738) 


0.016 


0 


0.703 


0.420 


Latin American 


allele frequency 
comparison 


A vs. G 


0.048 (0.003-0.826) 


0.036 


0 


1.000 


# 




additive model 


AG vs. GG 


0.058 (0.003-1 .034) 


0.053 


0 


1.000 


# 




additive model 


AA vs. GG 


0.291 (0.012-7.235) 


0.452 


0 


1.000 


# 




dominant model 


AA+AG vs. GG 


0.001 (0.001-0.008) 


0.001 


0 


1.000 


# 




recessive model 


AA vs. AG+GG 


0.314 (0.013-7.797) 


0.480 


0 


1.000 


# 



CI, confidence interval; 12, inconsistency index; OR, odds ratio; PEXG, pseudoexfoliation glaucoma; PEXS, pseudoexfoliation syndrome; SNP, single nucleotide 
polymorphism; vs., versus. Bold text: significant odds ratio (P<0.05) and significant between-study heterogeneity (P<0.1). 
# Publication bias could not be tested because a minimum of 3 studies were required. 
doi:1 0.1 371/joumal.pone.0090331 .t003 

P LOS ONE | www.plosone.org 5 March 2014 | Volume 9 | Issue 3 | e90331 



Association between L0XL1 and PEXS/PEXG 



Model 


Group by 


Study name 




Ethnicity 






Asian 


Rampnasad, 2008 




Asian 


Lee, 2009 




Asian 


Abu-Amero, 2010 




Asian 


Micheal. 2012 




Asian 


Park, 2013 


Fixed 


Asian 




Random 


Asian 






Caucasian 


Fan, 2008 




Caucasian 


Pasutto, 2008 




Caucasian 


MaluMewicz, 2011 




Caucasian 


MetaxaM, 2013 


Fixed 


Caucasian 




Random 


Caucasian 






Latin American 


Jaimes!2012 


Fixed 


Latin American 




Random 


Latin American 




Fixed 


Overall 




Random 


Overall 





Statistics for each study 



Odds 


Lower 


Upper 






ratio 


limit 


limit 


Z-Value 


p-Value 


0.135 


0.055 


0.328 


-4.413 


0.000 


0.003 


0.000 


0.024 


-5.579 


0.000 


0.025 


0.009 


0.068 


-7.251 


0.000 


0.022 


0.010 


0.051 


-8.999 


0.000 


0.005 


0.001 


0.021 


-7.159 


0.000 


0.029 


0.018 


0.047 


-14.598 


0.000 


0.019 


0.006 


0.064 


-6.502 


0.000 


0.009 


0.003 


0.027 


-8.449 


0.000 


0.026 


0.018 


0.038 


-19.327 


0.000 


0.004 


0.000 


0.080 


-3.660 


0.000 


0.309 


0.132 


0.722 


-2.710 


0.007 


0.033 


0.024 


0.046 


-20.746 


0.000 


0.030 


0.006 


0.149 


-4.294 


0.000 


0.000 


0.000 


0.008 


-5.250 


0.000 


0.000 


0.000 


0.008 


-5.250 


0.000 


0.000 


0.000 


0.008 


-5.250 


0.000 


0.031 


0.024 


0.040 


-25.741 


0.000 


0.015 


0.006 


0.038 


-9.040 


0.000 



Odds ratio and 95% O 



■ m 






<. 






c 






c 






< 






i 






i 






i 






i 






(. 








■ 




i 






i 






< 












< 






< 






< 







0.1 0.2 0.5 1 2 5 10 

AA+AG GG 



Figure 3. Forest Plots for associations between rs3825942 A/G (AA+AG vs. GG) and PEXS/PEXG stratified by ethnicity. The first author 
of the study and year of publication are shown for each citation. The squares and horizontal lines correspond to the study specific odds ratio (OR) and 
95% confidence interval (CI); the diamonds represent the pooled OR and 95% CI. 
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Table 4. Summary of pooled odds ratios for the association of rs21 65241 and PEXS/PEXG in the meta-analysis. 
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CI, confidence interval; 12, inconsistency index; OR, odds ratio; PEXG, pseudoexfoliation glaucoma; PEXS, pseudoexfoliation syndrome; SNP, single nucleotide 
polymorphism; vs., versus. Bold text: significant odds ratio (P<0.05) and significant between-study heterogeneity (P<0.1). 

# Publication bias could not be tested because a minimum of 3 studies were required. 

* Only allele frequency data were extracted from Ozaki et al. (2008). 

^The data for minor homozygosity+heterozygosity versus major homozygosity were extracted from the study of Park et al. 
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tibility to PEXS/PEXG was significant in all genetic models (Table 4). Furthermore, the subgroup analysis showed that all 
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Figure 4. Forest Plots for associations between rs21 65241 C/T (CC+CT vs. TT) and PEXS/PEXG stratified by ethnicity. The first author 
of the study and year of publication are shown for each citation. The squares and horizontal lines correspond to the study specific odds ratio (OR) and 
95% confidence interval (CI); the diamonds represent the pooled OR and 95% CI. 
doi:10.1371/journal.pone.0090331.g004 



ethnicities had significant associations for all of the genetic models 
(Table 4) and reduced susceptibility to PEXS/PEXG (Figure 4). 

Between-study heterogeneity analysis 

The QHest suggested significant between-study heterogeneity 
for several of the pooled models for each of the SNPs (Tables 2, 3, 
and 4). To examine the possible sources of heterogeneity, the 
studies were stratified by ethnicity, but the inconsistency index did 
not substantially decrease. We therefore performed a meta- 
regression analysis to assess the influence of the study character- 
istics on the effect estimates. 

The male proportion and mean age of the subjects did not 
significandy affect the influences of rs 1048661 and rs2 165241 on 
PEXS/PEXG susceptibility. In contrast, for the effects of 
rs3825942 on PEXS/PEXG risk in Caucasians, a significant 
influence was detected for both the male proportion (Figure 5A; 
slope: 0.272; 95% CI: 0.167-0.376; P = 0.0001) and mean age 
(Figure 5B; slope: 0.796; 95% CI: 0.375-1.217; P = 0.0002) of the 
subjects. 

Sensitivity analyses 

Although the genotype distribution in four of the studies did not 
follow Hardy- Weinberg equilibrium [16,27,28,29], the corre- 
sponding pooled ORs were not materially altered with or without 
including these studies (data not shown). In addition, the sensitivity 
analyses indicated that no single study had undue influence on the 
pooled OR results. 

Publication bias 

A Begg funnel plot analysis and Egger tests were used to access 
the publication bias of the literature. The shapes of the funnel plots 
of dominant models for the SNPs did not reveal evidence of 
obvious asymmetry (Figure SI). The Egger test results suggested 
that no publication bias was found in any of the comparison 
models (Tables 2, 3, and 4). 



Discussion 

PEXS/PEXG is a disorder characterized by the accumulation 
of abnormal fibrillar deposits in the anterior segment of the eye. 
LOXL1, which serves both as a cross-linking enzyme and an 
element of the scaffold to ensure spatially defined deposition of 
elastin and collagen substrates [7], was recently identified by 
genetic linkage studies as associated with a susceptibility to PEXS/ 
PEXG [14,15,31]. A large number of high-frequency SNPs have 
been identified for this gene [11,26]. After screening these SNPs 
with the selection criteria, the rsl048661 T/G, rs3825942 A/G, 
and rs2 165241 C/T polymorphisms were chosen to examine for 
their association with PEXS/PEXG susceptibility in this meta- 
analysis. 

Our results indicated that rsl048661 TT carriers had 92.1% 
and 40.4% less risk of developing PEXS/PEXG than did the 
controls in Caucasian and Asian populations, respectively, but had 
no influence on the susceptibility in the Latin American 
population. Carriers of rs3825942 AA or rs2165241 CC also 
had significantly less risk of developing PEXS/PEXG than did the 
non-carriers. Despite the ethnic heterogeneity in the LOXL1 
genotypes and the consequent variable susceptibility to PEXS/ 
PEXG, the rsl048661, rs3825942 and rs2165241 SNPs may 
provide a powerful diagnostic tool to identify subjects who are 
more likely to develop PEXS/PEXG. 

The analyses for heterozygosity of the variants are probably due 
to chance; the studies with small sample sizes for the minor 
homozygous alleles in PEXS/PEXG subjects would have insuf- 
ficient statistical power to detect slight effects, or they may have 
generated a fluctuated risk estimate. Given this situation, our 
evaluation of the effects associated with heterozygosity for the 
variants in our analysis should be interpreted with caution. 

In meta-analysis studies, heterogeneity could potentially restrict 
the interpretation of the pooled estimates, and ethnicity could play 
a role in determining the heterogeneity among studies. The 
different allele frequencies among the different ethnicities were a 
strong cause of the heterogeneity, leading us to do a subgroup 
analysis stratified by ethnicity; however, this resulted in no 
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Figure 5. Association between study population characteristics and effect of rs3825942 on PEXS/PEXG susceptibility in Caucasian 
patients. 

doi:10.1371/journal.pone.0090331.g005 



substantial decrease in the heterogeneity. Although applying a 
random effect model allowed us to estimate the effects of the 
different studies, we also performed a meta-regression to assess the 
influence of the population characteristics on the effect estimates. 
The results showed that the male proportion and mean age of the 
PEXS/PEXG patients positively correlated with the effect of SNP 



rs3825942 on the PEXS/PEXG susceptibility in the Caucasian 
population [32]. Although our study cannot explain how LOXL1 
interacts with gender and mean age, the correlations between 
them in the Caucasian population could be the primary cause of 
the heterogeneity. Considering the LOXL1 SNP genotypes 
together with these clinical predictors (gender and mean age) 
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may allow for greater accuracy in predicting the probability of 
developing PEXS/PEXG. 

To our knowledge, this study is the first meta-analysis to assess 
the association of the LOXL1 polymorphisms with PEXS/PEXG. 
This statistical method increased the power to detect and quantify 
an effect, and it provided a control for population differences that 
could lead to spurious associations if there are differences in gene 
frequency among the groups. Furthermore, this method allowed us 
to confirm the reliability and stability of the meta-analysis by 
performing publication bias and sensitivity analyses. 

Some limitations of this study should be taken into consider- 
ation. First, the study populations were primarily Caucasian and 
Asian. The subgroup meta-analysis for ethnicity had little or no 
information for other ethnic groups. Thus, strengthening the 
statistical power will require more data from other ethnic groups. 
Second, although we were able to discern a significant association 
between the LOXL1 polymorphisms and susceptibility to PEXS/ 
PEXG in the Caucasian and Asian populations, the sample size 
after pooling the existing studies was still relatively small. Third, 
the lack of available data prevented an adjustment for subgroup 
factors such as age, gender, and other variables that can interact 
with genetic factors to influence the marginal association estimates 
between the SNPs and PEXS/PEXG susceptibility. Should such 
data become available, a more precise analysis allowing for the 
adjustment of other covariates such as age, family history, 
environmental factors, and lifestyle would be feasible. 
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